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The objectives of the present work are twofold. The first is to address and resolve some of 
^SJ ' the differences present in independent, chiral-effective-field-theory (xEFT) derivations up to one 

T-H , loop, recently appeared in the literature, of the nuclear charge and current operators. The second 

■ objective is to provide a complete set of xEFT predictions for the structure functions and tensor 
04 ■ polarization of the deuteron, for the charge and magnetic form factors of ^He and "^H, and for the 

^ \ charge and magnetic radii of these few-nucleon systems. The calculations use wave functions derived 

^ . from high-order chiral two- and three-nucleon potentials and Monte Carlo methods to evaluate the 

■ relevant matrix elements. Predictions based on conventional potentials in combination with xEFT 
I charge and current operators are also presented. There is excellent agreement between theory and 

lO . experiment for all these observables for momentum transfers up to q < 2.0-2.5 fm~^; for a subset of 

■ them, this agreement extends to momentum transfers as high as g ~ 5-6 fm~^. A complete analysis 
' ^ '' of the results is provided. 

> . 

I . PACS numbers: 12.39.Fe, 13.40.-f 

Ph ■ I. INTRODUCTION 

T-H ■ Nuclear electromagnetic charge and current operators in chiral effective field theory (^EFT) up to one loop were 
^ ] derived originally by Park et al. [l] in covariant perturbation theory. Recently, two independent derivations, based on 
'jf^ ■ time-ordered perturbation theory (TOPT), have appeared in the literature, one by some of the present authors 0, Q 
] and the other by KoUing et al. [1, 0] . The expressions in Refs. 0, g and g for the tvifo-pion-exchange charge and 
I ■ current operators are in agreement with each other. Differences between the expressions reported in Refs. [H, Is'l and 
those in Ref. Q are found in some of the loop corrections to the one-pion-exchange (OPE) and short-range currents as 

■ well as the minimal currents originating from four-nucleon contact interactions involving two gradients of the nucleon 

■ fields. The differences in the loop corrections have their origin in the different implementations of TOPT adopted in 
. Refs. [1, [1] and Ref. Q, and relate to the treatment of reducible diagrams. One of the objectives of the present work 

• • ■ is to resolve some of these differences. This is addressed in Sec. Ulland Appendices \K\ and [Bl 
. ^ '. The other objective is to provide predictions for the charge and magnetic radii and form factors of the deuteron and 

■ trinucleons (^He and ^H), by utilizing two- and three-nucleon potentials derived either in xEFT or in the conventional 
5^ ■ framework, in combination with the charge and current operators obtained here. The methods used to carry out the 
5^ . calculations are discussed in Sec. IIIIl and a detailed analysis of the results is presented in Sec. IIVI This last section 

is organized into three subsections: the first illustrates the different strategies adopted for the determination of the 
low-energy constants (EEC's) that characterize the current operator up to one loop (no unknown EEC's enter the 
one-loop charge operator); the second and third report results, respectively, for the A(ci) and B{q) structure functions 
and tensor polarization T2Q{q) of the deuteron, and for the charge and magnetic form factors of ■^He and ^H, as well 
as results for the charge and magnetic radii of these few-nucleon systems. The conclusions are summarized in Sec. IVl 
while details on the evaluation of the loop integrals entering the charge operator are relegated in Appendix [C] 

There have been earlier xEFT studies of the deuteron electromagnetic structure in Refs.[^-Q and, most recently, 
in Ref. Q — this latter work has focused on the B{q) structure function. To the best of our knowledge, however, the 
one-loop xEFT predictions reported here for the ^He and '^H elastic form factors are new. 



II. NUCLEAR CHARGE AND CURRENT OPERATORS UP TO ONE LOOP 



The two-nucleon current (j) and charge (p) operators have been derived in xEFT up to one loop (to order eQ) in 
Refs. and respectively. In the following, we denote the momentum due to the external electromagnetic field 



2 



with q, and define 



k, = p^-p, , K, = (p^+p0/2, (2.1) 
k=(ki-k2)/2, K = Ki+K2, (2.2) 



where pi (p^) is the initial (final) momentum of nucleon i. We further define 

+1 +1 



n— — 2 n— — 3 

where the superscript n in j^"' and p^""^ specifies the order eQ" in the power counting. The lowest-order (LO) 
contributions j^"^-* and p*-^'^-' consist of the single-nucleon current and charge operators, respectively: 

j^"^' - — ^ [2eAr.i(g^)Ki + i^Ar,i(g2)cri X q] 

X5(p^-P2) + 1^2 , (2.4) 

and 

= eejv,i(g')<5(P2-P2) + 1^2 , (2.5) 

where rriN is the nucleon mass, q = k^ with z = 1 or 2 (the (5-functions enforcing overall momentum conservation 
q = ki have been dropped for simplicity here and in the following). 



eN.iil 



2, _ G|(g2) + 6-^(^2), 
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^^N,M') = ^'Miq') + GW)n^. (2.6) 

S /V S /V 

and and G^,| denote the isoscalar/isovector combinations of the proton and neutron electric (E) and magnetic 

(M) form factors, normalized as G|(0) = G')^(O) = 1, Gf^(O) = 0.880 ^jv, and GIj{0) = 4.706 /xjv in units of the 
nuclear magneton /z^v- The counting eQ~^ i'^Q~^) of the leading-order current (charge) operator results from the 
product of a factor eQ (eQ°) due to the coupling of the external electromagnetic field to the individual nucleons, 
and the factor from the momentum i5-function entering this type of disconnected contributions. Of course, this 
counting ignores the fact that the nucleon form factors themselves also have a power series expansion in Q. Here, 
they are taken from fits to elastic electron scattering data off the proton and deuteron p^— specifically, the Hohler 
parametrization Tl'| — rather than derived consistently in chiral perturbation theory (xPT) The calculations of 
the A = 2 and 3 nuclei elastic form factors that follow are carried out in the Breit frame, in which the electron- 
energy transfer vanishes. Hence, the hadronic electromagnetic form factors are evaluated at four-momentum transfer 

At order n = — 1 (NLO) there is a one-pion exchange (OPE) contribution to the current operator which reads 



ki 



X ■ 



(T2 • k; 



2 



K2 



1^2, (2.7) 



where we have defined uj'^ — k + rn^, TOtt being the pion mass. However, there are no n = —2 contributions to the 
charge operator. The presence of the isovector electric form factor G\ in j^^^-* follows from the continuity equation 



(2.8) 



where [. . . , . . .] denotes the commutator, is the charge operator given in Eq. ()2.5p . and v'"^^ is the static OPE 

potential 

4"nk) = -||r,.r2^i:^. (2.9) 
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FIG. 1: Diagrams illustrating the two-body charge operators at order n = or eQ". Nucleons, pions and photons are denoted 
by solid, dashed, and wavy lines, respectively. The solid circle in panel (a) is associated with a ^nN vertex of order eQ. Only 
one among the possible time orderings is shown. 



The l.h.s. of Eq. (|2.8p is of order Q^, the same as the r.h.s. since the commutator brings in an additional factor due 
to the implicit momentum integrations. It should be emphasized that the continuity equation requires that the same 
form factor be used to describe the electromagnetic structure of the hadrons in the longitudinal part of the current 
operator and in the charge operator. However, it places no restrictions on the electromagnetic form factors which 
may be used in the transverse parts of the current. Ignoring this ambiguity, the choice made here (G)^) satisfies the 
"minimal" requirement of current conservation [l3l |. 

Relativistic corrections to the leading order one-body current and charge operators enter, respectively, at n = 
and n = — 1 (both denoted as N2L0), and are given by 



i(o) = 



8m% 



ejv,i(g') 2 (K^ + qyA) (2Ki 



-l-i CTi X q) + Ki • q (q + 2 i CTi X Ki) 



"'"-AT 

X (4 (Ti X Ki - i q) - (2 i Ki - (Ti x q) g^/2 
+2 (Ki X q) cri • Ki] Sip'^ - p2) + 1 ^ 2 , 

+2 i q • (Ti X Ki) 5{p'^ - P2) + 1 ^ 2 , 
while the n = (N3L0) OPE two-body charge operators, illustrated in Fig.[l] read 



-,(-1) 



-,(0) 



9l 



(Ti • q cr2 • k2 



[G%{q^)T,-T,+GW)T2.] 

1^2 , 



g\ (Ti • k2 <r2 • k2 



(1-^) 



X [G-Uq^)Ti-T2+Gl{q^)T2,,] q-k2 
+ 2 I Gl{q^) (Xi X k2 • [ (1 - i^) Ki 

+ (1+;^) K2] 



1^2, 



.(0) 



9l 



(Ti • ki 0-2 • k2 



1 ^ 2 



(2.10) 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



The operator of panel (a) is due to a ^ttN vertex of order e Q originating from the interaction Hamiltonian 



2 tjinFt^ 



j dxA^V- (Vyl°) (r •7r-|-7r^)iV 



derived first by Phillips In the context of meson-exchange phenomenology, an operator of precisely this form 
results from considering the low-energy limit of the relativistic Born diagrams associated with virtual pion photo- 



4 



production amplitudes, see the review paper 13] and references therein. From this perspective, it appears reasonable 
to include the nucleon form factors and G]^ in Eq. (|2.12p . 

The operator of panel (b) depends on the off-energy-shell extrapolation, specified by the parameter adopted for 
the non-static corrections of order to the OPE potential [isj . 



(2.15) 



As shown in Ref. [l^ (and within the present approach in Ref. Q), different off-shell prescriptions for i;(^)(i/) and 
p''^^{v) are unitarily equivalent: 



(-3) 



(2.16) 



where the hermitian operator U{v) admits the expansion 

U{v) = u'-°\iy) + U'^^\v) + . . . , (2.17) 

and U''^''{v) and U^^\h') (see below) have been constructed, respectively, in Rcfs. and fsj (in this last paper, 
Eqs. (28) and (55), which give equivalent momentum-space expressions for U^^\v)^ contain a typographical error: the 
imaginary unit on the l.h.s. should be removed). Phenomenological potentials, such as the Argonne wis (AV18) [l^, 
and xEFT potentials, such as those recently derived by Entem and Machleidt [l3|, make the choice v = 1/2 in 
Eq. (|2.15l) . i.e., ignore non static corrections to the OPE potential. 

The operator of panel (c) , containing the ^tttt vertex, is obtained by expanding the energy denominators as Q 



1 



Ei — Ej 



1 



1 



Ei — El 



(2.15 



where Ei denotes NN (or NN^) intermediate energies and the pion energy (or energies, as the case may be), and 
by noting that the leading (static) corrections vanish, when summed over the possible six time orderings. However, the 
terms proportional to the ratio {Ei — i?/)/w^, which is of order Q, lead to the non-static operator given in Eq. (|2.14p . 
It is multiplied by the pion form factor G^((7^), which we parametrize in vector-meson dominance and consistently 
with experimental data at low momentum transfers as 



1 



(2.19) 



where nip is the p-meson mass. 



A. Current operators at order n = 1 (eQ) 

The currents at order eQ (N3L0) are illustrated diagrammatically in Fig. [51 and consist of: (i) terms generated by 
minimal substitution in the four-nucleon contact interactions involving two gradients of the nucleon fields as well as 
by non-minimal couplings to the electromagnetic field; (ii) OPE terms induced by jnN interactions beyond leading 
order; and (iii) one-loop two-pion-exchange (TPE) terms. We discuss them below. 

The contact minimal and non minimal currents, denoted by the subscripts "min" and "nm" respectively, are written 

as 



5(1) 



I e 



Ja.nm 



Gliq^) {Tl X T2), [(G2 + 3 G4 + Gt) ki 

+(G2 - G4 - G7)ki (Tl • 0-2 

+G7(Ti • (ki - k2) 0-2 - ^ eAr,i(g^) G5 

x((Ti + (T2) X ki + 1 ^ 2 , 



X q+ 1 ^ 2 . 



(2.20) 



(2.21) 
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FIG. 2: Diagrams illustrating the two-body current operators at order n — 1 oi eQ. Nucleons, pions and photons are denoted 
by solid, dashed, and wavy lines, respectively. The solid circle in panel (b) is associated with a jttN vertex of order eQ^. Only 
one among the possible time orderings is shown. 

The expression above for j^^nin Fierz-transformed version of the current given in Eq. (3.11) of Ref. 0, see 

App. 1^ for a derivation. We note that the first three terms in Eq. (|2.20p agree with the first line of Eq. (5.3) of 
Kolling et al. j5|, while the term proportional to C5 differs by the isoscalar piece, which, however, can be absorbed 
in a redefinition of C[^. The low-energy constants (EEC's) Ci, . . . ,Ct, which also enter the two-nucleon contact 
potential, have been constrained by fitting np and pp elastic scattering data and the deuteron binding energy. We 
take their values from the Machleidt and Entem 2011 review paper The EEC's CJg and C(g (and dg, d'g, and 
(^21 below) are determined by fitting measured photo-nuclear observables of the A — 2 and 3 systems, as discussed 
in Sec. IIVI Finally, we observe that there is no a priori justification for the use of G^/G^ (or Gff/G\j) in the non- 
minimal contact current, and these form factors are included in order to provide a reasonable fall-off with increasing 
for the strength of this current. 

The isovector (IV) OPE current at N3L0 is given by 



= (1) 

Jb.IV 



'd'2i{Ti X T2)z (Ti X k2 



X q+ 1 ^ 2 



(2.22) 



and depends on the two (unknown) EEC's dg and (^21. They can be related [2] to the A^-A transition axial coupling 
constant and magnetic moment (denoted as ^^na) in a resonance saturation picture, which justifies the use of the 
jNA electromagnetic form factor for this term. It is parametrized as 



GjNAiq'^ 



fJ"yNA 



(l + 9VAi,i)yi + gVAl2 



(2.23) 



where h-jna is taken as 3 /i at from an analysis of jN data in the A- resonance region [l8[ . This analysis also gives 
Aa,i=0.84 GeV and Aa, 2=1-2 GeV. The isoscalar (IS) piece of the OPE current depends on the LEC d'g mentioned 
earlier. 



.(1) . 9a ^ , 2\ 
Jb,IS = * e «9 G^^p{q ) T-i • T2 



T2 • k2 



k2 X q + 1 ^ 2 



(2.24) 



k-2 



and, again in a resonance saturation picture, reduces to the well known ^irp current [2]. Accordingly, we have 
accounted for the q^ fall-off of the electromagnetic vertex by including a ^irp form factor, which in vector-meson 
dominance is parametrized as 



GfTrpiq'^) — 



1 -I- (p- jm^ 



(2.25) 



rrii^ is the w-meson mass. We can now clarify the differences in these tree-level currents as reported here and in 
Ref. We first note that the relations between the primed d[ and di in Ref. ^] should have read: dg — — 8dg, 
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dg = — Sdg, d.2i = 2(^21 — ^22- The term proportional to (^22 originates from the Lagrangian 

term oc d22 - |^ ^22 ^^^'^ [^''(t x tt), F^,] TV . 
Kolhng et al. f§\ integrate it by parts to obtain 



(2.26) 



2e 



term oc ^22 = — d22 



+N^S''{TXTv),F^4d''N) , 



(2.27) 



while the authors of Ref. use the equations of motion for the electromagnetic field tensor at leading order, d^F^^^ — 0, 
to express it as 



term oc ^22 = |^ ^22 F^,^, (r x d^Tv), N. 



(2.28) 



These two different treatments lead to the ^22 current as given in Ref. [2;] and [5| . They differ by a term proportional to 
(cTi X q) X q, which does not contribute to the magnetic moment (Ml) operator fi — —{i/2) x j|g_g 0- Similarly, 

the term proportional to f5(q) in Eq. (4.28) of Ref. [5] does not give any contribution to /i, since f5{q) <x. for small q. 
The term proportional to f%{q) is included in Eq. (|2.7p provided g\ — > g\{l — dig m^/gA). The value adopted here 
for gA is obtained from two-nucleon scattering data (Sec. IIV[) . Therefore, for processes induced by Ml transitions, 
such as the n d and n ^He radiative captures at thermal neutron energies studied in Ref. (loj or the magnetic scattering 
under consideration in this work, the differences above are irrelevant. 
The one-loop TPE currents, diagrams (c)~(k) of Fig. [21 are written as 



J loop 



I e Gl{q^) (n X X2), Vfe Fi(fc) + I e Gl{q^) T2,, 
k (Ti • k 



Fo{k)(T^-F2{k)^ 



fc2 



X q+ 1 



(2.29) 



where the functions Fi{k) are 



Foik) 



Fi(fc) 



F2{k) 



9\ 



^9a 
fc2 + 4 m?. 



G{k) 



4(l + gi)m2 16 glmi 



k^ +4 ml 
^ G(fc) 



(fc2 + 4 774)2 



96 7r2Fi 



2-2<?l 



Amiil 



-k\l + 10gl-23g\)- 



1%g\mi 



9\ 



2-65! 



fc2 + 4 m 



4 + k"^ 

G{k) 



''9a ml 



45i 



4(l + 35i)m2 16gim4 



fc2 4- 4 TO?. 



(fc2 + 4to2)2 



(2.30) 



(2.31) 



(2.32) 



and the loop function G(/c) is defined as 



G(fc) 



k'^ J\ ml + k'^ + k 
— m ■ 



fc2 



(2.33) 



The expression above results from expanding j{oQp(q, k) in a power series in q as j{oop(q, k) = JioopC^'^) — *q x 
/i.^''''(k) + where /x^^^ is the magnetic dipole operator, and j[QQp(0, k), which corresponds to the first term in 
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Eq. (|2.29l) . satisfies current conservation with the TPE potential v^J(k) (of order Q^), since 



4'J(k-q/2),eAra 



-1^2 



-ie(xi X X2)2q- VfcFi(fc) + 1^2 



(2.34) 



to leading order in q. In fact, the current jioop(0,k) is proportional to the electric dipole operator, and does not 
contribute to elastic electromagnetic transitions, such as those of interest here. 

Finally, we note that a more careful analysis, detailed in Appendix|Bl of the loop short-range currents corresponding 
to diagrams (h)-(k) in Fig.[2]shows that they vanish, in contrast to that which was reported in Ref. and in agreement 
with the result of Ref. H . 



B. Charge operators at order n = 1 (e Q) 

The two-body charge operators at one loop (N4L0) are illustrated in Fig. [31 and have been derived in Ref. [3]. The 
contributions from diagrams of type (a)-(b) and (g)-(h) vanish, and after carrying out the loop integrations (discussed 
in App. ICj). those from diagrams of type (c)-(f) and (i)-(j) read: 



•,(1) 



1 o2 nl/2 



4L(x,fc2) 



1 2 



pi'' 



J-i^Gliq') / d.. 



L{x,k2) 

1 ^ 2 



4L{x, fci) 



L{x,ki) 

[4t2^2 + (ti X T2)J 



(2.35) 



(2.36) 



-24L(x,fc2) ' 



L{x,k2) L^{x,k2) 



Pi 



(1) 



[4ti^^ - i/ (ti X T2)^] 



((72 X k2) ■ ((Ti X k2 

L(x, ^2) 



1 ^ 2 



(2.37) 









dx X dy 




-15\ix,y) + 


'0 J-1/2 







3A-(B + C) + (A + B)-(A + C) 



X{x,y) 

-{(Ti X A) • {(T2 X A) - (cTi X A) • {(T2 X C) - (cTi X B) • ((T2 X A) + (cTi X B) • (0-2 X C) 



1 



(A • B)(A • C) + (Ti • (A X B) (T2 • (A X C 

1 



X^{x,y) 

-B • (0-2 X A) (A + B) • ((T2 X C 



1 



A2(a;,y) 



X{x,y) 
A-Bcr2-(Ax C) 



(ti X T2)z 

hl^2 



3 0-2 • (A X C) 



(2.38) 



Pi 



(1) 



1 5^ 

e--4-CTG)^((?^)ri,2cri • (T2 + 1 ^ 2 



-e-^CTG^{q^)T,,J dx 
-^Fi Jo 
1/4 -a;2 



3L2(x,q)-TO2 



L{x,q) 



-CTi ■ (T2 



L{x,q) 



CTi • q 0-2 • q 



1 ^ 2 



(2.39) 



(2.40) 



where we have defined 



L\x,p) = {l/A-x^)^^ + ml 



(2.41) 
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FIG. 3: Diagrams illustrating two-body charge operators entering at order n — 1 or eQ. Nucleons, pions, and photons are 
denoted by the solid, dashed, and wavy lines, respectively. Only one among the possible time orderings is shown. 



X'^{x,y) ^xq^/A- [x y q - (1 - x) k] % (1 - x)k^ + ml , 



(2.42) 



A = -x{yq + k) , (2.43) 
B = {l-2xy)q/2 + (1- x) k , (2.44) 
C = -(l + 2xy)q/2+(l-a;) k . (2.45) 

It is easily verified that the charge operators (c)+(d), (c) + (f), and (i)+(j) vanish at q = 0. Finally, we note that the 
form of the operator (e) depends on the ofF-the-energy-shell prescription adopted for the non-static corrections to the 
TPE potential. As in the OPE case, however, these different forms for the TPE non-static potential and accompanying 
charge operator are unitarily equivalent 3], in particular Pc^\iy) = Pc^\o) + [p'-^'^-' , iU^^\h')'\. In closing, we note 
that the re-analysis, outlined in Appendix |B1 of the loop corrections to the short-range charge operators illustrated 
in panels (g)-(j) has led to expressions which are different from those reported originally in Ref. They also differ 
from those in Ref. 

III. CALCULATION 

The deuteron charge (Gc), magnetic (Gm), and quadrupole (Gq) form factors are obtained from [S^l 

d-M) , (3.1) 

d;0)] , (3.2) 



(3.3) 

where | d; M) is the deuteron state with spin projection — M , p and jy denote, respectively, the charge operator 
and y component of the current operator, the momentum transfer q is taken along the z-axis (the spin quantization 
axis), and rj = [q/2mdY (™d is the deuteron mass). They are normalized as 

Gc(0) = 1 , GAf(O) - (md/mA,) pa , Gq(0) = , (3.4) 

where pd and Qd are the deuteron magnetic moment (in units oi pn) and quadrupole moment, respectively. Expres- 
sions relating the form factors to the measured structure functions A and -B, and tensor polarization T2Q are given 



Gc{q) - o E ^d-M\p{qi) I 

M=±lfl 

Gniq) = -^lia[{d;l\jy{qi)\ 
V ^ V 

Ggiq) = ^[{d;0\p{qz)\d;0) 
-{d;l\ p{qz) I d;l)] , 



9 



in Ref. [20j. The calculations are carried out in momentum space [20| with techniques similar to those described in 
some detail below for the trinucleons. 

The charge and magnetic form factors of the trinucleons are derived from 

Fc{q) = |(+|p(gz)|+) , (3.5) 

FM{q) = -^Im[(-h,(gz)|+)] , (3.6) 

with the normalizations 

Fc(0) = l, FM(0)=/i, (3.7) 

where /i is the magnetic moment (in units of /iat). Here | ±) represent either the '^He state or •^H state in spin 
projections Jz = ±1/2. In momentum space, the one-body electromagnetic operators in Sec. HIl have the generic form 

Oib(q) = - ^(^™) ^(^") Oib(k,,K0 , (3.8) 

cyclic l,m,n 

and their matrix elements can be written as 



(Oib(q)) = X! / V'lf'(Pi +q/2,Pm,p^ 



cyclic/, m,n 

xOib(q,Pi)^M(pz -q/2, ) , (3.9) 

where we have defined 

'^'^^ and S{---) = (2^)^^(---) • (3-10) 



(2 7r)3 



For an assigned configuration (p;, p,„, Pn), the wave functions are expanded on a basis of 8 x 3 spin-isospin states for 
the three nucleons as 



24 



l/'(p;,p,„,p„) = ^ Va(P!,Pm,Pn) | a) , (3.11) 

a=l 

where the components ipa are complex functions and the basis states (for '^H, for example) | a) =\ {p t)i, [n t)2, {n t)3)i 
I (n t)i, (p 1)2, {n 1)3), and so on. The spin-isospin algebra for the overlaps 

24 

^to^^ ^ raOabA , (3.12) 
a, 6=1 

is carried out with the techniques developed in Ref. [21|]. Monte Carlo (MC) methods are used to evaluate the 
integrations in Eq. p.9p by sampling momenta from a (normalized) probability density | ipM{Pi,Pm,Pn) P according 
to the Metropolis algorithm. 

The two-body operators in Sec. HIl have the momentum-space representation 

02b (q) = ^(K;™ - q) ^(k„) 

cyclic l,m.n 

x02b(Ki„/2 + ki™,Ki™/2- W) , (3.13) 

where the momenta Kjm = k; -|- km and k/„i = (k; — km)/2. These operators have power law behavior at large 
momenta, and need to be regularized. This is accomplished by introducing a momentum cutoff function of the form 

CAihm) = e-(^--/^)' , (3.14) 
with the parameter A in the range (500-600) MeV (see discussion in Sec. lIVp. The matrix elements are expressed as 
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(02b(q)) = E / / Vl/'(Pi+q/4 + k,„,/2,p„ + q/4-k,„/2,p„)CA(fc,™) 

cyclic ;,m,n-^'^''" "'Pi.P^.Pr. 

x02b(q,k/„)V'A/(P/-q/4-k,„/2,p,„-q/4 + k,,„/2,p„) . (3.15) 

The spin-isospin algebra is handled as above, while the multidimensional integrations are efficiently done by a com- 
bination of MC and standard quadratures techniques. We write 

(02b(q))-/dk/ F{k,pi,pm,Pn)^^Y.WU^ ^ (3-16) 
J Jp,,p„,p„ ^^c^^^yv\c) 

where c denotes configurations (k,p/,Pm,p„) (total number Nc) sampled with the Metropolis algorithm from the 
probability density W{c) —\ V'j\/(P;, Pm, Pn) P/(47r), i.e., uniformly over the k directions. For each such configuration 
c, the function F is obtained by Gaussian integration over the magnitude kim (as well as the parameters x and y for 
the case of the charge operators at one loop) 



I r°° 

^(c) = Y] 77r\3 dfc,„fcf„ V Vl(---fc;mk...) 



cyclic l.m,n a, 6—1 

x02hMb{q,kl7nk)lpb{- ■ -kniii-- ■ ■) ■ (3-17) 

Convergence in these Gaussian integrations requires of the order of 20-30 points, in the case of kim distributed over 
a non-uniform grid up to 2 A or so, while Nc of the order of 100, 000 is sufficient to reduce the statistical errors in 
the MC integrations, which are of the order of a few % at the highest q values (and considerably smaller at lower 
q). These MC errors are further reduced by taking appropriate linear combinations of the matrix elements of the 
electromagnetic operators using different q directions and different spin projections for the initial and final states. The 
trinucleons wave functions are obtained with the hyperspherical harmonics (HH) expansion discussed in Refs. [22l-[2^. 
This method can be applied in either coordinate- or momentum-space. Below, we briefly review its momentum-space 
implementation. 

A. The hyperspherical harmonics method in momentum-space 

The trinucleon wave functions with total angular momentum JJz are written as 

\^''~')=J2'^^\<'') ' (3-18) 

where IV'^''^) is a suitable complete set of states, and /x is an index denoting the set of quantum numbers necessary to 
specify the basis elements (see below). By applying the Rayleigh-Ritz variational principle, the problem of determining 
and the ground-state energy Eq of the system is reduced to a generalized eigenvalue problem. 
In momentum space we define the Jacobi momenta as 



k2p = {Pj - P^) /V2 , kip = ^/2/3[pk - ip^ + Pj) /2] , (3.19) 

where denotes the momentum of nucleon i and p specifies a given permutation of the three nucleons, with p — 1 
corresponding to the ordering 1,2,3. We introduce a hyper-momentum K and a set of angular and hyper-angular 
variables as 

K - (kl^ + fc2p) , nf^ = [k2p, kip ; ^p] , (3.20) 
where tan^p = kip/k2p . In terms of these variables, the basis functions IV'^'^"') are defined as 

\^'/'^) = gGi{K)y{G}{n^'''>) , (3.21) 
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where y{G}{^'-^^) are written as 1251 



3 

E 

P=i L 



Y, 



LL 

[G] 



















SS^_ 



T2 



(3.22) 



and the sum is over the three even permutations. The spins (isospins) of nucleons i and j are coupled to S2 {T2), 
which is then coupled to the spin (isospin) of the third nucleon to give a state with total spin S (isospin TTz). The 
total orbital angular momentum L and total spin S are coupled to the total angular momentum JJj. The functions 



Yr^f with definite values of LLz are the hyperspherical- harmonics functions, and are written as [22 



'[G] 



LL 

[G] 



YtA^2p 



5^^. (kip) 

LL^ 



N[G] (cOS(/)p) 

(cos2(/)p) , 



(3.23) 



where Y^j(kip) and Yi^(\!i2p) are spherical harmonics, N^q^ is a normalization factor, and + (cos2(/)p) denotes 

the Jacobi polynomial of degree n. The grand angular quantum number G is defined as G = 2 rt+i?i+i?2. The subscripts 
{G} and [G] in Eqs. stand, respectively, for {G} = {£1, £2, i, S'2, T2, 5, T; n} and [G] = [^i,€2;ri], and /i 

in Eq. (I3.18P stands for ji = {G} I. Finally, the functions gGi{K) in Eq. (13.211) are defined as 



9Gi{K) = 



dpp^ Jg+2{Kp) flip) , 



(3.24) 



where Jg+2 [Kp] are Bessel functions and the functions /; [p) are related to Laguerre polynomials {ip) via 

(5), 



m = 7V^!/a + 5)! Lf\^p) e--^"'^ 



(3.25) 



The non-linear parameter 7 is variationally optimized. With this form of /; {p) , the corresponding functions gc i (K) 
can easily be calculated, and are explicitly given in Ref. [23| . The form adopted for gGi{K) is such that the momentum- 
space basis is simply the Fourier transform of the coordinate-space one |24| . 



IV. RESULTS 



This section consists of three subsections. In the first one, we discuss various strategies for the determination of the 
unknown LEC's dg, dg, d^i^ G(5, and G(g entering the current operator at N3L0. In contrast, the charge operator 
up to N4L0 only depends on the nucleon axial coupling constant gA, pion decay amplitude Pjr, and nucleon mass 
and magnetic moments. The values adopted in the present work for gA and Ft, are, respectively, 1.28 and 184.6 
MeV, which give a TriV coupling constant [g^^j^i^) of 13.6, as obtained in analyses of NN elastic scattering data at 
energies below the pion production threshold [26| . The two-body operators are regularized via the cutoff function in 
Eq. and A values of 500 MeV and 600 MeV are considered. 

In the second and third subsections we present results, respectively, for the deuteron A{q) and B{q) structure 
functions and tensor polarization T2a{q), and for the charge and magnetic form factors of '^H and '^He, along with 
results for the static properties of these few- nucleon systems including the deuteron quadrupole moment, the deuteron 
and trinucleons charge and magnetic radii and magnetic moments. The A = 2 calculations use either the Argonne 
?;i8 (AV18) [Hi or chiral potentials at order with cutoff set at 500 MeV (N3L0) or 600 MeV (N3L0*) [T?]. Of 
course, the A = 3 calculations also include three-nucleon potentials — the Urbana-IX model ^2^ in combination with 
the AV18, and the chiral N2L0 potential ^ in combination with either the N3L0 or N3L0*. The LEC's cd and ce 
(in standard notation) in the chiral three-nucleon potential have been constrained by reproducing the ^H/'^He binding 
energies and the tritium Gamow- Teller matrix element [2§| in each case. With the AV18/UIX Hamiltonian, the 
and ^He binding energies are found to be 8.487 MeV and 7.747 MeV, respectively. 

The calculations are carried out in configuration space in the first subsection, and in momentum space — with the 
methods outlined in Sec. IIIH — in the following two subsections. We have checked that the r- and p-space versions 
of the computer codes produce identical results up to to tiny differences due to numerics and to numerically non- 
equivalent implementations of the momentum cutoff function in Eq. (I3.14[) in these r- and p-space calculations. The 
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TABLE I: Adimensional values of the isoscalar LEC's corresponding to cutoffs A = 500 MeV and 600 MeV obtained for the 
N3LO/N2LO and N3LO*/N2LO* Hamiltonians; the values in parentheses are relative to the AV18/UIX Hamiltonian. 



A 


df 


df X 10 


500 
600 


4.072 (2.522) 
11.38 (5.238) 


2.190 (-1.731) 
3.231 (-2.033) 



TABLE II: Cumulative contributions to the deuteron and trinucleons isoscalar magnetic moments in units of /ijv, corresponding 
to cutoffs A = 500 MeV and 600 MeV obtained for the N3LO/N2LO and N3LO*/N2LO* Hamiltonians; the contributions in 
parentheses are relative to the AV18/UIX Hamiltonian. The experimental values for the deuteron and trinucleons isoscalar 
magnetic moments are 0.8574 /iiv and 0.4257 /i]v, respectively. 







MS 


A 


500 


600 


500 


600 


LO 
N2L0 
N3L0(min) 
N3LO(nm) 
N3L0(0PE) 


0.8543 (0.8472) 
0.8471 (0.8400) 
0.8725 (0.8739) 
0.8548 (0.8593) 
0.8574 (0.8574) 


0.8543 (0.8472) 
0.8474 (0.8400) 
0.8806 (0.8760) 
0.8538 (0.8626) 
0.8574 (0.8574) 


0.4222 (0.4104) 
0.4143 (0.4027) 
0.4501 (0.4455) 
0.4247 (0.4269) 
0.4257 (0.4257) 


0.4220 (0.4104) 
0.4155 (0.4027) 
0.4611 (0.4483) 
0.4235 (0.4313) 
0.4257 (0.4257) 



hadronic electromagnetic form factors entering the one- and two-body charge and current operators are those specified 
in Sec. |lll The matrix elements of these operators are evaluated in the Breit frame with Monte Carlo methods. The 
number of sampled configurations is of the order of 10^ for the deuteron and 10^ for the A — 3 systems. The statistical 
errors, which are not shown in the results that follow, are typically < 1% over the whole momentum-transfer range, 
and in fact much less than 1% for q < 2 fm~^. 

A. Determination of the LEC's 

As already remarked, the LEC's Ci, i = 1,...,7, in the minimal contact current, corresponding to A cutoffs of 
500 and 600 MeV, are taken from fits to NN scattering data [l7|. In reference to the LEC's entering the OPE and 
non-minimal contact currents at N3L0, it is convenient to introduce the adimensional set (in units of the cutoff 
A) as 

C[, = df/A^ , 4 = dl/A^ 

c[, = dY/A^ , 4 = , 4i = > 

(4.1) 

s V 

where the superscript or V on the d^ ' characterizes the isospin of the associated operator, i.e., whether it is isoscalar 
or isovector. The isoscalar df , listed in Table HI have been fixed by reproducing the experimental deuteron magnetic 
moment /i^ and isoscalar combination /is of the trinucleon magnetic moments. The LEC df multiplying the contact 
current is rather large, but not unreasonably large, while the LEC df is quite small. The cumulative contributions 
to fj,d and are reported in Table |TT1 The NLO and N3LO-loop magnetic moment operators are isovector, and 
therefore do not contribute to these isoscalar observables. At N3L0 the only non-vanishing contributions are those 
associated with the OPE and minimal (min) and non-minimal (nm) contact currents. Of course, the last row in 
Table HIl reproduces the experimental values for /i^ and /is. 

The isovector LEC d^ is taken as (i^/4 by assuming A dominance. The three different sets of remaining LEC's dY 
and d^ reported in Table Hill have been determined in the following way. In set I d^ and d^ have been constrained to 
reproduce the experimental values of the np radiative capture cross section anp at thermal neutron energies and the 
isovector combination fiy of the trinucleons magnetic moments. This procedure, however, leads to unreasonably large 
values for both LEC's, and is clearly unacceptable. In particular, it makes the contributions of the associated magnetic 
dipole operators unnaturally large, and, as shown in Table IIVI totally spoils the expected convergence pattern. This 
pathology is especially severe in the case of the AV18/UIX Hamiltonian model. 

In sets II and III d^ is assumed to be saturated by the A resonance, i.e. 
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TABLE III: Adimensional values of the isovector LEC's corresponding to cutoffs A — 500 MeV and 600 MeV obtained for the 
N3LO/N2LO and N3LO*/N2LO* Hamiltonians; the values in parentheses are relative to the AV18/UIX Hamiltonian. Note 
that da = in all cases; see text for further explanations. 



A 


dY{i) 




dYiii) 


^^(11) 


dr(III) 


dxm 


500 
600 


10.36 (45.10) 
41.84 (257.5) 


17.42 (35.57) 
33.14 (75.00) 


-13.30 (-9.339) 
-22.31 (-11.57) 


3.458 
4.980 


-7.981 (-5.187) 
-11.69 (-1.025) 


3.458 
4.980 



TABLE IV: Cumulative contributions to the np radiative capture cross section in mb and trinucleons isovector magnetic 
moment in units of /iiv, corresponding to cutoffs A = 500 MeV and 600 MeV obtained for the N3LO/N2LO and N3LO*/N2LO* 
Hamiltonians; the contributions in parentheses are relative to the AV18/UIX Hamiltonian. See text for further explanations. 
The experimental values for the np cross section and trinucleons isovector magnetic moment are (332.6±0.7) mb and —2.553 hn, 
respectively. 









A 


500 


600 


500 


600 


LO 
NLO 
N2LO 
N3L0(loop) 
N3L0(min) 


305.8 (304.6) 
320.6 (319.3) 

319.2 (317.7) 

321.3 (320.9) 
321.3 (320.9) 


304.6 (304.6) 
318.9 (320.9) 
317.6 (319.2) 
320.5 (322.4) 
320.5 (322.4) 


-2.193 (-2.159) 
-2.408 (-2.382) 
-2.384 (-2.359) 
-2.430 (-2.418) 
-2.413 (-2.406) 


-2.182 (-2.159) 
-2.392 (-2.413) 
-2.370 (-2.390) 
-2.432 (-2.448) 
-2.415 (-2.437) 


N3L0(nm, <-I) 
N3LO(OPE, dX-1) 


315.2 (287.4) 
332.6 (332.6) 


305.7 (242.7) 
332.6 (332.6) 


-2.297 (-1.782) 
-2.553 (-2.553) 


-2.142 (-0.9029) 
-2.553 (-2.553) 


N3LO(nm, df-II) 
N3L0(0PE, d^-II) 


329.1 (328.1) 
332.6 (332.6) 


328.5 (326.2) 

332.6 (332.6) 


-2.562 (-2.535) 
-2.612 (-2.610) 


-2.561 (-2.506) 
-2.622 (-2.616) 


N3L0(nm, ^-111) 
N3LO(OPE, d^-III) 


326.0 (324.9) 
329.4 (329.4) 


324.7 (322.7) 

328.8 (329.1) 


-2.502 (-2.478) 
-2.553 (-2.553) 


-2.491 (-2.443) 
-2.553 (-2.553) 



where ttia ~ rn^ — 294 MeV, hA/F.^ — f-^NAlmiiT, with f^NA/i^^) = 0.35 as obtained by equating the first-order 
expression of the A-decay width to the experimental value, and the transition magnetic moment h^na = 3 fiN — a 
similar strategy has been implemented in a number of calculations, based on the xEFT magnetic moment operator 
derived in Ref. [l|, of the np, nd, and n'^He radiative captures, and magnetic moments of ^ = 2 and 3 nuclei (soj . On 
the other hand, the LEG dY multiplying the contact current is fitted to reproduce either cr„p in set II or /iy in set III. 
Both alternatives still lead to somewhat large values for this LEG, but we find the degree of unnaturalness tolerable 
in this case. We observe that there are no three-body currents at N3L0 [l^, and therefore it is reasonable to fix the 
strength of this Ml operator by fitting a three-nucleon observable such as /iy. 

Gumulative contributions to (Jnp and /iy are listed in Table HVl At N3L0, we have identified separately those due only 
to loop currents labeled as N3LO(loop), and those from loop+minimal contact currents labeled as N3L0(min). The 
experimental values for CT„p and fiv are reproduced with set I, row labeled N3L0(0PE, d2-I), while only or 
are reproduced with set II or III, rows labeled N3L0(0PE, d^-II) or N3L0(0PE, rf^-III). Indeed, the N3L0(0PE, 
d^-II or III) results provide predictions for fiy or anp, respectively. These predictions are within 3% for fxy and 1% 
for (T„p of the experimental values, and exhibit a weak cutoff and Hamiltonian-model dependence. 

In Ref. the df'^ were determined using the same procedure adopted here for set I. However, the values reported 
in that work are drastically different from those obtained in the present one. These differences are due to several 
factors: i) in Ref. [l^ the Ml operator derived from Eq. (|2.29p included an isovector loop correction proportional 
to the leg's Cs and Ct, which turns out to vanish in a more careful analysis of the relevant diagrams (the loop 
short-range currents discussed in Appendix |B| ; ii) in Ref. [l^ the values for the LEGs Ci , . . . , C7 were taken from a 
chiral potential obtained at (NLO) Q rather than at (N3L0) [l3| as in the present case; iii) in Ref. [l^ the 
minimal contact current is the Fierz-transformed version of that given in Eq. (j2.20l) (see discussion in Appendix |^ . 
However, this Fierz equivalence is spoiled by the regularization procedure, i.e. by the inclusion of the same cutoff 
function Cf^{k) for both. Hence the contribution of this current in the present work is different from that obtained in 
Ref. [3. 

B. Static properties and form factors of the deuteron 
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TABLE V: Cumulative contributions to the deuteron root-mean-square charge radius and quadrupole moment corresponding 
to cutofTs A — 500 and 600 MeV obtained with the N3LO and N3LO* Hamiltonians; results in parentheses are relative to the 
AV18 Hamiltonian. The experimental values for rd and Qd are 1.9734(44) fm 31] and 0.2859(3) fm'^ 32], respectively. 





rd (fm) 


Qd (fm^) 


A 


500 


600 


500 


600 


LO 
N2LO 
N3LO(OPE) 
N3LO(i^ = 1/2) 


1.976 (1.969) 
1.976 (1.969) 
1.976 (1.969) 
1.976 (1.969) 


1.968 (1.969) 
1.968 (1.969) 
1.968 (1.969) 
1.968 (1.969) 


0.2750 (0.2697) 
0.2731 (0.2680) 
0.2863 (0.2818) 
0.2851 (0.2806) 


0.2711 (0.2697) 
0.2692 (0.2680) 
0.2831 (0.2814) 
0.2820 (0.2802) 



The deuteron root-mean-square charge radius and quadrupole moment, obtained with the chiral and AV18 potentials 
and cutoff parameters A = 500 MeV and 600 MeV, are listed in Table IVl We denote the leading order (n = — 3 in 
the notation of Sec. HT]) term of Eq. (|2.5I) with LO, the n = — 1 relativistic correction of Eq. (|2.11l) with N2L0, and 
the n = terms of Eqs. ([2T2|) and (|2T3l) ~ ([2T4l) with N3L0(0PE) and N3L0(i/), respectively. The remaining charge 
operators at N4L0 {n = 1), being isovector, do not contribute to these observables (and corresponding form factors). 
The N3LO/N3LO* and AV18 potentials neglect retardation corrections in their OPE component, which corresponds 
to setting v = 1/2 in Eq. (12.15^ . Note that the isoscalar piece of the N3L0(j^) charge operator scales as 1 — i^, 
and contributes less than 0.5% of the LO result for = 1/2. The N2L0 and N3L0 corrections to rd, which is well 
reproduced by theory, are negligible. The chiral potential predictions for Qd are within 1% of the experimental value, 
while the AV18 ones underestimate it by about 2%. Variation of the cutoff in the (500-600) MeV range leads to about 
1% (negligible) changes in the N3LO/N3LO* (AVIS) results. The LO and N2L0 charge operators do not include the 
cutoff function and the AV18 results are independent of A. This is not the case for the results corresponding to the 
N3L0 and N3L0* potentials because of their intrinsic A dependence. 



The deuteron A{q) structure function and tensor polarization T'2o(<z), obtained at LO and by including corrections 
up to N3L0 in the charge operator, are compared to data in Fig. SI top panels. In this figure (as well as in those that 
follow) the momentum-transfer range goes up to g = 7.5 fm~^, much beyond the ~ 3-4 m-^ upper limit, where one 
would naively expect this comparison to be meaningful, given that the present theory retains up to TPE mechanisms. 
On the other hand, we note that the next (non- vanishing) isoscalar contributions only enter at N5L0 {n = 2) Q , and 
are therefore suppressed by two powers of Q relative to those at N3L0. 

The A{q) structure function is well reproduced by theory up to g ~ 3 fm^^. At higher momentum transfers, the 
N3L0 results based on the AV18 tend to overestimate the data — a feature also seen in the conventional approach of 
Ref. 11^ — while those based on the chiral potentials still provide a good fit to the data. The cutoff dependence is 
weak at low g, but becomes more pronounced as q increases. 

Similar considerations hold for the T2o(g) observable, although in this case the N3L0 results derived from the 
chiral potentials overpredict the data for g > 3 fm~^^^hile those from the AV18 fit reasonably well the data up to 
q ~ 4.5 fm~^. In contrast, the conventional approach [20[ (also based on the AV18, of course) reproduces very well the 
measured T20 over the whole g-range. The OPE charge operator in that work has the same structure as the present 
N3L0(0PE) one, but includes a much harder cutoff than adopted here. Furthermore, the calculation of Ref. [20] also 
retains short-range (isoscalar) mechanisms associated with p-meson exchange and ^irp transition, which in ^EFT are 
presumably subsumed in contact operators at N5L0 [sj. We note that in both ^(g) and T2o(g) a small magnetic 
contribution, discussed separately below, is accounted for (the electron scattering angle in T20 is set at 70°). 

The charge and quadrupole form factors extracted from the unpolarized and tensor polarized deuteron data are 
compared to results obtained in LO and by including corrections up to N3L0 in Fig. SI bottom panels. The Gc{q) and 
Gqlq) form factors calculated with deuteron wave functions from the chiral potentials are in qualitative agreement 
with predictions obtained by Phillips Q at the same chiral order (although the potentials used in that study 
are from Ref. [l^ rather than from Ref. [l3| as in the present work). The spread in the N3L0 results due to cutoff 
variations observed here is similar to that reported in Ref. @ for both Gc{q) and GQ{q). However, the central values 
for these observables in the momentum-transfer region g > 2.5 fm~^ reported in that work appear to underestimate 
the data appreciably. This is not the case here, particularly for GQ(g), for which the N3L0 predictions provide an 
excellent fit to the measured values (up to g ~ 6 fm~^). These differences likely arise from differences in the deuteron 
wave functions obtained in Ref. flT^ and Ref. |59| (see Fig. 16 in the 2011 review paper T?] for a comparison). Indeed, 
for these same reasons, the AVIS results are in better agreement with data for Gc{q) in the diffraction region than 
the N3LO/N3LO*, while the reverse is true for Gqiq) at g > 3 fm"^ The AVIS deutero n wave function, particularly 




FIG. 4: (Color online). The deuteron A{q) structure function and tensor polarization T2o{q) (top panels), and charge and 
quadrupole form factors Gc(g) and GqIq) (bottom panels), obtained at leading order (LO) and with inclusion of charge 
operators up to N3LO (TOT), is compared with experimental data from Refs. |33454|. Predictions corresponding to = 1/2 
and cutoffs A in the range 500-600 MeV are displayed by the bands. 

its D-wave component, is markedly different from the N3LO/N3LO* (see again Fig. 16 in the Machleidt and Entem 
review [TTj). 

The individual contributions corresponding to A = 500 MeV and the N3L0 potential are listed in Tables IVll and Ivnl 
for q values in the range (0.0-5.1) fm~^. The N2L0 (N3L0) charge operators are proportional to ^/mj^ (I/totv), 
and therefore vanish in the static limit. The N3L0(0PE) correction is the leading one for q > 1.5 fm^^, and is 
responsible for shifting the zero in the LO Gc{q) to lower q. However, this correction interferes constructively with 
the LO contribution in the case of GQ{q). The i^-dependent retardation correction N3L0(i^) is found to be negligible, 
which allows one to conclude that violations of the unitary equivalence between the OPE potential and associated 
charge operator is of little numerical import (for v — 0-1). 

The deuteron magnetic moment is one of the two observables utilized to fix the LEC's entering the isoscalar current 
operators at N3L0, denoted as N3L0(nm) and N3L0(0PE) in Sec. HV^ and Table HD The structure function B{q) 
and magnetic form factor Gm{<i), obtained with the AV18 and chiral potentials, and currents at LO and by including 
corrections up to N3L0, are compared to data in Fig. [51 There is generally good agreement between theory and 
experiment for q values up to ~ 2 fm^^. At higher g's, the results corresponding to the chiral (AV18) potential 
under-predict (over-predict) the data significantly when the current includes up to N3L0 corrections. In particular, 
the diffraction seen in the data at g ~ 6.5 fm^^ is absent in the AV18 calculations, and is shifted to lower q values in 
the N3LO/N3LO* ones. There are large differences between the N3LO/N3LO* and AVIS results with the LO current, 
which simply reflect differences in the S- and D-wave components of the deuteron wave functions corresponding to 
these potentials. The cutoff dependence is large for the chiral potentials, while it remains quite modest for the AVIS 
over the whole momentum transfer range. This is consistent with the rather different sensitivity of the LEC's df 
and df to variations of A in the (500-600) MeV range obtained with either the chiral potential or AVIS, see Table ID 
There is a mismatch in the chiral counting between the potentials of Ref. [l3| at order and the present current 
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TABLE VI: Individual contributions to the monopole form factor Gc{q) corresponding to cutoff A = 500 MeV for the N3L0 
Hamiltonian; v = 1/2 and (— x) stands for 10~^. 



q(fm-i) 


LO 


N2L0 


N3LO(OPE) 


N3LO( 


^) 


0.0 


1.00 




0.00 




0.00 


0.00 




0.3 


0.945 




-0.340( 


-3) 


-0.211(-3) 


-0.600( 


-6) 


0.6 


0.792 




-0.113( 


-2) 


-0.799(-3) 


-0.500( 


-6) 


0.9 


0.614 




-0.196( 


-2) 


-0.165(-2) 


0.530( 


-5) 


1.2 


0.452 




-0.255( 


-2) 


-0.260(-2) 


0.218( 


-4) 


1.5 


0.321 




-0.280( 


-2) 


-0.349(-2) 


0.515( 


-4) 


1.8 


0.220 




-0.273( 


-2) 


-0.422(-2) 


0.932( 


-4) 


2.1 


0.146 




-0.242( 


-2) 


-0.470(-2) 


0.143( 


-3) 


2.4 


0.920(- 


-1) 


-0.197( 


-2) 


-0.493(-2) 


0.195( 


-3) 


2.7 


0.547(- 


-1) 


-0.145( 


-2) 


-0.493(-2) 


0.245( 


-3) 


3.0 


0.295(- 


-1) 


-0.923( 


-3) 


-0.474(-2) 


0.287( 


-3) 


3.3 


0.131(- 


-1) 


-0.448( 


-3) 


-0.441(-2) 


0.319( 


-3) 


3.6 


0.295(- 


-2) 


-0.518( 


-4) 


-0.400(-2) 


0.339( 


-3) 


3.9 


-0.278(- 


-2) 


0.250 


-3) 


-0.356(-2) 


0.348( 


-3) 


4.2 


-0.556(- 


-2) 


0.453 


-3) 


-0.311(-2) 


0.346( 


-3) 


4.5 


-0.645(- 


-2) 


0.566 


-3) 


-0.269(-2) 


0.335( 


-3) 


4.8 


-0.621(- 


-2) 


0.602 


-3) 


-0.230(-2) 


0.318( 


-3) 


5.1 


-0.539(- 


-2) 


0.579 


-3) 


-0.196(-2) 


0.297( 


-3) 



TABLE VII: Same as in Table IVTl but for the quadrupole form factor GQ{q), normalized at g = as in Eq. H3.4[) . 



q(fm-^) 


LO 


N2LO 


N3L0(0PE) 


N3L0( 




0.0 


24.8 


-0.172 




1.20 


-0.108 




0.3 


23.5 


-0.176 




1.19 


-0.986( 


-1) 


0.6 


19.7 


-0.182 




1.13 


-0.968( 


-1) 


0.9 


15.3 


-0.184 




1.04 


-0.936( 


-1) 


1.2 


11.4 


-0.179 




0.930 


-0.890( 


-1) 


1.5 


8.29 


-0.167 




0.810 


-0.829( 


-1) 


1.8 


5.92 


-0.150 




0.690 


-0.757( 


-1) 


2.1 


4.18 


-0.131 




0.576 


-0.677( 


-1) 


2.4 


2.92 


-0.111 




0.473 


-0.596( 


-1) 


2.7 


2.03 


-0.917(- 


-1) 


0.383 


-0.515( 


-1) 


3.0 


1.39 


-0.740(- 


-1) 


0.307 


-0.440( 


-1) 


3.3 


0.947 


-0.584(- 


-1) 


0.244 


-0.371( 


-1) 


3.6 


0.635 


-0.449(- 


-1) 


0.193 


-0.310( 


-1) 


3.9 


0.418 


-0.335(- 


-1) 


0.152 


-0.257( 


-1) 


4.2 


0.268 


-0.242(- 


-1) 


0.119 


-0.212( 


-1) 


4.5 


0.167 


-0.168(- 


'1) 


0.933(-l) 


-0.174( 


-1) 


4.8 


0.997(-l) 


-0.111(- 


-1) 


0.731(-1) 


-0.142( 


-1) 


5.1 


0.568(-l) 


-0.686(- 


-2) 


0.574(-l) 


-0.116( 


-1) 



at order eQ. This becomes obvious when considering current conservation, which for these potentials would require 
accounting for terms up to order eQ^ in the current, well beyond available derivations 0,0,0] at this time. 



The AV18 results obtained here for B{q) are similar to those reported in the conventional framework of Ref. [l^] 
(see curve labeled IA+p7r7-NR in Fig. 5). In that work, the current included the standard impulse-approximation 
(lA) term — the LO current in ^EFT — and the two-body term from pnj transitions. The size, and in fact sign, of 
the pnj contribution were found to depend on whether the current was derived by retaining the fully relativistic (R) 
structure of the associated Feynman amplitude, or only the leading-order term in its non-relativistic (NR) expansion — 
in this latter case, it is essentially the N3L0(0PE) current of Eq. (|2.24p . Indeed, the pn^ contribution had the same 
(opposite) sign as the lA when it was evaluated with the NR (R) current, and the IA-|-p7r7(NR) results overestimated 
the data by an amount similar to that shown in Fig. [5] 




FIG. 5: (Color online). The deuteron B{q) structure function (top panel) and magnetic form factor GM{q) (bottom panel), 
obtained at leading order (LO) and with inclusion of current operators up to N3L0 (TOT) , is compared with the experimental 
data from Refs. [3a . [39l liol . [ssI - IstI ]. Predictions corresponding to cutoffs A in the range 500-600 MeV are displayed by the 
bands. 



TABLE VIII: Individual contributions to the magnetic form factor Gm{<i), normalized at 5 = as in Eq. H3.4|l . corresponding 
to cutoff A = 500 MeV for the N3LO Hamiltonian; {-x) stands for lO""". 



q(fm-^) 


LO 


N2L0 


N3LO(min) 


N3LO( 


nm) 


N3LO(OPE) 


0.0 


1.71 




-0.144( 


-1) 


0.508 


-1) 


-0.353 


-1) 


0.483( 


-2) 


0.3 


1.62 




-0.150( 


-1) 


0.508 


-1) 


-0.353 


-1) 


0.480( 


"2) 


0.6 


1.36 




-0.159( 


-1) 


0.495 


-1) 


-0.343 


-1) 


0.473( 


-2) 


0.9 


1.07 




-0.166( 


-1) 


0.474 


-1) 


-0.329 


-1) 


0.461( 


-2) 


1.2 


0.807 




-0.166( 


-1) 


0.447 


-1) 


-0.310 


-1) 


0.444( 


-2) 


1.5 


0.590 




-0.159( 


-1) 


0.416 


-1) 


-0.289 


-1) 


0.423( 


-2) 


1.8 


0.422 




-0.144( 


-1) 


0.382 


-1) 


-0.265 


-1) 


0.399( 


-2) 


2.1 


0.296 




-0.125( 


-1) 


0.347 


-1) 


-0.241 


-1) 


0.372( 


-2) 


2.4 


0.202 




-0.104( 


-1) 


0.312 


-1) 


-0.217 


-1) 


0.344( 


-2) 


2.7 


0.134 




-0.818( 


-2) 


0.279 


-1) 


-0.194 


-1) 


0.315( 


-2) 


3.0 


0.860(- 


-1) 


-0.608( 


-2) 


0.248 


-1) 


-0.172 


-1) 


0.286( 


-2) 


3.3 


0.522(- 


-1) 


-0.418( 


-2) 


0.219 


-1) 


-0.152 


-1) 


0.258( 


-2) 


3.6 


0.291(- 


-1) 


-0.254( 


-2) 


0.193 


-1) 


-0.134 


-1) 


0.231( 


"2) 


3.9 


0.138(- 


-1) 


-0.120( 


-2) 


0.169 


-1) 


-0.117 


-1) 


0.206( 


-2) 


4.2 


0.415(- 


-2) 


-0.159( 


-3) 


0.148 


-1) 


-0.103 


-1) 


0.183( 


-2) 


4.5 


-0.152(- 


-2) 


0.581 


-3) 


0.130 


-1) 


-0.899 


-2) 


0.162( 


-2) 


4.8 


-0.444(- 


^2) 


0.105 


-2) 


0.113 


-1) 


-0.785 


-2) 


0.143( 


-2) 


5.1 


-0.554(- 


-2) 


0.128 


-2) 


0.988 


-2) 


-0.686 


-2) 


0.126( 


-2) 




FIG. 6: (Color online). The ^He and charge form factors (top panels), and their isoscalar and isovector combinations 
(bottom panels), obtained at leading order (LO) and with inclusion of charge operators up to N4L0 (TOT), is compared with 
experimental data [53|. Predictions corresponding to v = 1/2 and cutoffs A in the range (500-600) MeV are displayed by the 
bands. 

Recently, a calculation of the deuteron magnetic structure, based on the same ^EFT utilized here, has appeared in 
the literature Q . It uses chiral potentials at order derived in Ref. fE^ , and a different strategy from that adopted 
here for constraining the two LEC's in the isoscalar N3L0 current. One of them is still fixed by reproducing /i^ ; the 
other, however, is determined by a fit to Bi^q) data up to g ~ 2 fm^^. Predictions for this observable in g = (2-4) 
fm~^ seem to overestimate the data at the highest q values (q > 3.5 fm~^), but display much less cutoff dependence 
than obtained here. This is clearly due to the different way in which the LEC's are constrained in the two calculations. 

Finally, in Table IVllIl we list the individual contributions to Guisi) obtained with the N3L0 potential and cutoff 
A = 500 MeV. The notation is as follows: LO is the leading-order (e Q^^) current of Eq. (|2.4p : N2L0 is the relativistic 
correction of order n = (e Q") in Eq. (piB : N3L0( min), N3L0(nm), and N3L0(0PE) are the corrections of order 
71 = 1 (eQ) in Eqs. ([2211), and dOll), respectively. The N3L0(min) and N3L0(nm) contributions from 

the minimal and non-minimal contact currents cancel to a large extent, and their combined effect is comparable to 
the N3L0(0PE) contribution. This interplay among different corrections, however, depends strongly on A and the 
Hamiltonian model considered. 



C. Static properties and form factors of the trinucleons 



The notation for the various components of the charge operator is the same as given at the beginning of Sec. IIVB[ 
except that now the one-loop (isovector) corrections at N4L0 contribute too, since the ^He and nuclei have 



19 



TABLE IX: Cumulative contributions in fm to the '^He and '^H root-mean-square charge radii corresponding to = 1/2 and 
cutofTs A = 500 MeV and 600 MeV, obtained with the N3LO/N2LO and N3L07N2LO* Hamiltonians; resuhs in parentheses 
are relative to the AV18/UIX Hamiltonian. The experimental values for the "^He and "^H charge radii are [sj (1.959 ± 0.030) 
fm and (1.755 ± 0.086) fm, respectively. 





^He 




A 


500 


600 


500 


600 


LO 
N2LO 
N3LO 
N4LO 


1.966 (1.950) 
1.966 (1.950) 
1.966 (1.950) 
1.966 (1.950) 


1.958 (1.950) 
1.958 (1.950) 
1.958 (1.950) 
1.958 (1.950) 


1.762 (1.743) 
1.762 (1.743) 
1.762 (1.743) 
1.762 (1.743) 


1.750 (1.743) 
1.750 (1.743) 
1.750 (1.743) 
1.750 (1.743) 




FIG. 7: The '^H charge form factor obtained with the N3LO/N2LO Hamiltonian, cutoff A = 500 MeV, and charge operators 
up to N4LO corresponding to = 0, 1/2, and 1. 



predominantly total isospin T = 1/2. As a matter of fact, the hyperspherical harmonics wave functions utilized to 
represent their ground states also include small T = 3/2 admixtures due to isospin-symmetry breaking terms induced 
by the electromagnetic and strong interactions. 

There are no unknown LEC's entering the charge operator up to N4L0, and the predicted root-mean-square 
charge radii of ^He and "^H, obtained with the N3LO/N2LO and AV18/UIX combinations of two- and three- nucleon 
potentials and cutoffs in the (500-600) MeV range, are listed in Table llXl Corrections at N2L0, N3L0, and N4L0 
are negligible — the corresponding operators vanish at g = 0. The spread between the N3LO/N2LO (A = 500 MeV) 
and N3L07N2LO* (A = 600 MeV) resuhs at LO is about 0.5%, which is much smaller, particularly for ^H, than 
the experimental error. The predicted radii for both Hamiltonian models are within 0.5% of the current experimental 
central values. 

The calculated charge form factors of ^He and '^H, and their isoscalar and isovector combinations F^{q) and {q), 
normalized, respectively, to 3/2 and 1/2 at g = 0, are compared to data in Fig. [SI The agreement between theory and 
experiment is excellent for q < 2.5 fm~^. At larger values of the momentum transfer, there is a significant sensitivity 
to cutoff variations in the results obtained with the chiral potentials. This cutoff dependence is large at LO and is 
reduced, at least in ^He, when corrections up to N4L0 are included. These corrections have opposite sign than the 
LO, and tend to shift the zeros in the form factors to lower momentum transfers, bringing theory closer to experiment 
in the diffraction region. 

As already remarked, the chiral (and conventional) two-nucleon potentials utilized in the present study ignore 
retardation corrections in their OPE and TPE components, which corresponds to the choice v = 1/2 in the non-static 
pieces of the corresponding potentials and accompanying charge operators in Eqs. (|2.13p and (|2.37p [3]. Figure [7] is 
meant to illustrate how inconsistencies between the potential and charge operator impact predictions for the form 
factor, by presenting results obtained with the N3LO/N2LO Hamiltonian (v — 1/2), cutoff A — 500 MeV, and N3L0 
and N4L0 corrections with = and 1 in the charge operator. Their effect is negligible. 
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FIG. 8: Cumulative contributions to the '^He and ''H charge form factors, obtained with the N3LO/N2LO Hamiltonian, cutoff 
A — 500 MeV, and u — 1/2, from the components of the charge operator order by order. 

In Fig. [3 we show cumulatively the LO, N2L0, N3L0, and N4L0 contributions to the charge form factors of '^He 
and ^H. The N2L0 are smallest, while the N3L0 and N4L0 turn out to be comparable. This is illustrated explicitly 
in Tables 1x1 and IXll where we list the individual contributions of the various terms entering at each order. In these 
tables, we denote with N3L0(0PE) and N3L0(z^) the operators in Eqs. (|2?T2|) and ([2J3l) - ([2T4l) . respectively; with 
N4L0(cd), N4L0(ef; v), N4L0(ij) those in Eqs. (jOHl) - (12:361) . Eqs. (|237l) - (12:381) . and Eqs. (l2:39| - (12:401) . Among 
the corrections at N3L0 the OPE term — column N3L0(0PE) — illustrated by panel (a) in Fig. [T] is dominant, while 
among those at N4L0 the TPE terms — columns N4L0(ef; ly) and N4L0(ij) — illustrated by panels (e)-(f) and (i)-(j) 
in Fig. |3| are dominant. The N3L0(0PE) and N4L0(ef; and N4L0(ij) contributions are of similar magnitude, 
indeed there is no hint of suppression in going from N3L0 to N4L0, as one would have naively expected on the basis 
of power counting. 

The '^He contributions in Table |X| have been divided by the number of protons Z — 2 in order to have the form 
factor normalized to one a.t q — 0. The N4L0 charge operators are isovector and, if '^He and '^H were pure T = 1/2 
states, then 2xN4LO('^He)=-N4LO(^H). That this equality is not exactly satisfied reflects the fact that the present 
•^He and wave functions are not simply the charge mirror of each other — that is, (J^^ Ti^x) pHe) 7^ pH), where Ti^x 
is the x-component of nucleon i isospin operator. 



Moving on to the magnetic structure of the trinucleons, we note that the isoscalar combination of ^He and 
magnetic moments is used to fix one of the two (isoscalar) EEC's entering the current at N3L0. Both the isovector 
combination /iy and the np radiative capture cross section (T„p are used to fix the isovector EEC's in set I of the 
N3L0 currents, while in sets II and III one of these EEC's is fixed by A dominance, and the other is determined by 
reproducing (T„p (mv) in set II (III), see Tables HIl and lIVl By construction, then, the ^He and magnetic moments 
are exactly reproduced in sets I and III, while in set II they are calculated to be, respectively, -2.186 (-2.196) /ijv and 
3.038 (3.048) with the N3LO/N2LO (N3LO*/N2LO*) Hamihonian and A = 500 (600) MeV, and similar resuhs 
with the AV18/UIX Hamiltonian. These should be compared to the experimental values of -2.127 /xjv and 2.979 /i^r. 
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TABLE X: Individual contributions to the '^He charge form factor, obtained with the N3LO/N2LO Hamiltonian, cutoff A — 500 
MeV, and u = 1/2; (-x) stands for lO""". 



q (fm-i) 


LO 


N2L0 


N3LO(OPE) 


N3LO(!/) 


N4LO(cd) 


N4LO(ef; 


N4L0(ij) 


0.0 


1.00 




0.00 




0.00 




0.00 




0.00 


0.00 


0.00 




0.2 


0.983 




-0.561(- 


-3) 


-0.152( 


-3) 


0.100(- 


-4) 


-0.477(-5) 


0.904(-4) 


0.301( 


-4) 


0.6 


0.807 




-0.406(- 


-2) 


-0.125( 


-2) 


0.870(- 


-4) 


0.108(-3) 


0.785(-4) 


0.346( 


-3) 


1.0 


0.562 




-0.760(- 


-2) 


-0.289( 


-2) 


0.205(- 


-3) 


0.296(-3) 


0.1 26 (-4) 


0.835( 


-3) 


1.4 


0.342 




-0.876(- 


-2) 


-0.434( 


-2) 


0.324(- 


-3) 


0.503(-3) 


-0.134(-3) 


0.132( 


-2) 


1.8 


0.185 




-0.758(- 


-2) 


-0.513( 


-2) 


0.409(- 


-3) 


0.680(-3) 


-0.347(-3) 


0.167( 


-2) 


2.2 


0.892(- 


-1) 


-0.525(- 


-2) 


-0.516( 


-2) 


0.447(- 


-3) 


0.796(-3) 


-0.585(-3) 


0.184( 


-2) 


2.6 


0.366(- 


-1) 


-0.288(- 


-2) 


-0.463( 


-2) 


0.433(- 


-3) 


0.841(-3) 


-0.796(-3) 


0.184( 


-2) 


3.0 


0.111(- 


-1) 


-0.107(- 


-2) 


-0.380( 


-2) 


0.385(- 


-3) 


0.821(-3) 


-0.938(-3) 


0.170( 


-2) 


3.4 


0.623(- 


-3) 


-0.969(- 


-5) 


-0.292( 


-2) 


0.315(- 


'3) 


0.750(-3) 


-0.989(-3) 


0.148( 


-2) 


3.8 


-0.258(-2) 


0.488(- 


-3) 


-0.213( 


-2) 


0.241(- 


-3) 


0.644(-3) 


-0.949(-3) 


0.122( 


-2) 


4.2 


-0.276(-2) 


0.577(- 


-3) 


-0.148( 


-2) 


0.174(- 


-3) 


0.522(-3) 


-0.831 (-3) 


0.957( 


-3) 


4.6 


-0.200(-2) 


0.476(- 


-3) 


-0.994( 


-3) 


0.117(- 


-3) 


0.395(-3) 


-0.660(-3) 


0.702( 


-3) 


5.0 


-0.119(-2) 


0.318(- 


-3) 


-0.642( 


-3) 


0.730(- 


-4) 


0.274(-3) 


-0.469(-3) 


0.476( 


-3) 



TABLE XI: Same as in Table H but for ^H. 



q (fm-i) 


LO 


N2LO 


N3L0(0PE) 


N3L0( 




N4LO(cd) 


N4L0(ef; u) 


N4LO(ij) 


0.0 


1.00 




0.00 




0.00 




0.00 




0.00 


0.00 


0.00 


0.2 


0.991 




0.647(- 


-3) 


-0.185( 


-3) 


0.110(- 


-4) 


0.103(-4) 


-0.190(-3) 


-0.633(-4) 


0.6 


0.844 




0.462(- 


-2) 


-0.152( 


-2) 


0.880(- 


-4) 


-0.229(-3) 


-0.136(-3) 


-0.729(-3) 


1.0 


0.621 




0.844(- 


-2) 


-0.356( 


-2) 


0.245(- 


-3) 


-0.628(-3) 


0.459(-4) 


-0.176(-2) 


1.4 


0.402 




0.940(- 


-2) 


-0.543( 


-2) 


0.462(- 


-3) 


-0.107(-2) 


0.396(-3) 


-0.278(-2) 


1.8 


0.231 




0.784(- 


-2) 


-0.653( 


-2) 


0.686(- 


-3) 


-0.144(-2) 


0.874(-3) 


-0.352(-2) 


2.2 


0.118 




0.520(- 


-2) 


-0.671( 


-2) 


0.858(- 


-3) 


-0.169(-2) 


0.138(-2) 


-0.389(-2) 


2.6 


0.517(- 


-1) 


0.270(- 


-2) 


-0.615( 


-2) 


0.944(- 


-3) 


-0.179(-2) 


0.182(-2) 


-0.389(-2) 


3.0 


0.174(- 


-1) 


0.903(- 


-3) 


-0.519( 


-2) 


0.940(- 


-3) 


-0.174(-2) 


0.210(-2) 


-0.361(-2) 


3.4 


0.204(- 


-2) 


-0.972( 


-4) 


-0.411( 


-2) 


0.867(- 


-3) 


-0.159(-2) 


0.218(-2) 


-0.314(-2) 


3.8 


-0.328(- 


-2) 


-0.483( 


-3) 


-0.310( 


-2) 


0.615(- 


-3) 


-0.136(-2) 


0.206(-2) 


-0.259(-2) 


4.2 


-0.403(- 


-2) 


-0.508( 


-3) 


-0.225( 


-2) 


0.543(- 


-3) 


-0.110(-2) 


0.178(-2) 


-0.201(-2) 


4.6 


-0.313(- 


-2) 


-0.383( 


-3) 


-0.159( 


-2) 


0.483(- 


-3) 


-0.825(-3) 


0.139(-2) 


-0.147(-2) 


5.0 


-0.196(- 


-2) 


-0.232( 


-3) 


-0.108( 


-2) 


0.364(- 


-3) 


-0.565(-3) 


0.942(-3) 


-0.985(-3) 



The ^He and '^H magnetic radii corresponding to sets I-III are given in Table IXlTl The predicted values are consistent 
with experiment, although the measurements have rather large errors (10% for '^H). Their spread as A varies in the 
(500-600) MeV range is at the 1% level or less. A recent quantum Monte Carlo study iGOll. using wave functions 
derived from conventional two- and three nucleon potentials (the AV18 and Illinois 7 model [6l|) and set III of xEFT 
currents, has led to predictions for magnetic moments and transitions in nuclei with A < 9 in excellent agreement 
with the measured values. Therefore in the following, unless stated otherwise, we adopt set III of isovector LEC's. 
We disregard set I for the reasons already explained in Sec. IIV Al 



TABLE XII: Cumulative contributions in fm to the ^He and root-mean-square magnetic radii corresponding to cutoffs 
A = 500 MeV and 600 MeV, obtained with the N3LO/N2LO and N3LO*/N2LO* Hamiltonians; results in parentheses are 
relative to the AV18/UIX Hamiltonian. Predictions corresponding to sets I, II, and II of isovector LEC's dY and dX in 
Table [Till are listed. The experimental values for the ^He and "^H magnetic radii are 57] (1.965±0.153) fm and (1.840±0.181) fm, 
respectively. 





^He 


«H 


A 


500 


600 


500 


600 


LO 
NLO 
N2LO 


2.098 (2.092) 
1.990 (1.981) 
1.998 (1.992) 


2.090 (2.092) 
1.983 (1.974) 
1.989 (1.984) 


1.924 (1.918) 
1.854 (1.847) 
1.865 (1.859) 


1.914 (1.918) 
1.845 (1.841) 
1.855 (1.854) 


N3L0(I) 


1.924 (1.931) 


1.910 (1.972) 


1.808 (1.800) 


1.796 (1.819) 


N3LO(II) 


1.901 (1.890) 


1.883 (1.896) 


1.789 (1.774) 


1.773 (1.778) 


N3L0(III) 


1.927 (1.915) 


1.913 (1.924) 


1.808 (1.792) 


1.794 (1.797) 




FIG. 9: (Color online). The "^He and "^H magnetic form factors (top panels), and their isoscalar and isovector combinations 
(bottom panels), obtained at leading order (LO) and with inclusion of current operators up to N3L0 (TOT) corresponding 
to the LEC's df and df in Table U and to set III of isovector LEC's dY and d^ in Table Hill are compared with experimental 
data [5^. Predictions relative to cutoffs A in the range (500-600) MeV are displayed by the bands. 

The magnetic form factors of ■^He and '^H and their isoscalar and isovector combinations F^j{q) and F]^j{q), nor- 
malized respectively as fxs and fiy at q = 0, at LO and with inclusion of corrections up to N3L0 in the current, 
are displayed in Fig. |9l As is well known from studies based on the conventional meson-exchange framework (see 
the review (isj ] and references therein), two-body currents are crucial for "filling in" the zeros obtained in the LO 
calculation due to the interference between the S- and D-state components in the ground states of these nuclei. For 
g < 2 fm~^ there is excellent agreement between the present ^EFT predictions and experiment. However, as the 
momentum transfer increases, even after making allowance for the significant cutoff dependence, theory tends to 
underestimate the data, in particular it predicts the zeros in both form factors occurring at significantly lower values 
of q than observed. Thus, the first diffraction region remains problematic for the present theory, confirming earlier 
conclusions derived from studies in the conventional framework (62l . l63j . 

Figure [TUl illustrates the sensitivity of the N3L0 predictions on the different ways in which the isovector LEC's are 
constrained in sets I, II, III. The set I results are strongly at variance with data. Set II leads to two-body current 
contributions larger than in set III, and consequently, in contrast to set III, the corresponding form factors reproduce 
the data in the diffraction region. However, the cutoff variation of the results is considerably larger than for set HI, 
as reflected in the change of the LEG dY for A = 500-600 MeV in Table HlTl Furthermore, set II overestimates nv by 
about 3%. 



Figure [TT] exhibits cumulatively the LO, NLO, N2L0, and N3L0 contributions to the "^He and '^H magnetic form 
factors, obtained with the N3LO/N2LO Hamiltonian and cutoff A = 500 MeV. Tables IXlIIl and IXTVl list the individual 
components of these contributions at selected values of q. The notation is as follows: with LO we denote the one-body 



23 




q(fm"') 



FIG. 10: (Color online). The ^He and magnetic form factors, obtained at leading order (LO) and with inclusion of current 
operators up to N3L0 (TOT) corresponding to sets I, II, and III of isovector LEC's dX and dX in Table Hill is compared with 
experimental data [S^. Predictions, relative to the N3LO/N2LO Hamiltonian and corresponding to cutoffs A = 500-600 MeV, 
are displayed by the bands. 



TABLE XIII: Individual contributions to the ''He magnetic form factor, obtained with the N3LO/N2LO Hamiltonian, cutoff 
A = 500 MeV, and set HI of isovector LEC's; {-x) stands for 10"^. 



q(fm-i) 


LO 


NLO 


N2LO 


N3L0 (L) 


N3L0 (min) 


N3LO(nm) 


N3L0(0PE) 


0.2 


-1.72 




-0.196 




0.162(- 


-1) 


-0.385( 


-1) 


0.511( 


-1) 


-0.109 




-0.488(-l) 


0.6 


-1.37 




-0.184 




0.160(- 


-1) 


-0.358( 


-1) 


0.478( 


-1) 


-0.101 




-0.414(-1) 


1.0 


-0.898 




-0.160 




0.142(- 


-1) 


-0.310( 


-1) 


0.421( 


-1) 


-0.877(- 


-1) 


-0.293(-l) 


1.4 


-0.495 




-0.128 




0.105(- 


-1) 


-0.253( 


-1) 


0.352( 


-1) 


-0.716(- 


-1) 


-0.164(-1) 


1.8 


-0.228 




-0.934(- 


-1) 


0.595(- 


-2) 


-0.196( 


-1) 


0.280( 


-1) 


-0.554(- 


-1) 


-0.584(-2) 


2.2 


-0.794( 


-1) 


-0.632(- 


-1) 


0.189(- 


-2) 


-0.145( 


-1) 


0.214( 


-1) 


-0.410(- 


-1) 


0.106(-2) 


2.6 


-0.964( 


-2) 


-0.402(- 


-1) 


-0.912(- 


-3) 


-0.103( 


-1) 


0.159( 


-1) 


-0.292(- 


-1) 


0.448(-2) 


3.0 


0.158(- 


-1) 


-0.241(- 


-1) 


-0.234(- 


-2) 


-0.709( 


-2) 


0.114( 


-1) 


-0.201(- 


-1) 


0.542(-2) 


3.4 


0.199(- 


-1) 


-0.138(- 


-1) 


-0.266(- 


-2) 


-0.473( 


-2) 


0.796( 


-2) 


-0.134(- 


-1) 


0.496(-2) 


3.8 


0.159(- 


-1) 


-0.756(- 


-2) 


-0.231(- 


-2) 


-0.306( 


-2) 


0.543( 


-2) 


-0.869(- 


-2) 


0.392(-2) 


4.2 


0.103(- 


-1) 


-0.396(- 


-2) 


-0.170(- 


-2) 


-0.191( 


-2) 


0.360( 


-2) 


-0.543(- 


-2) 


0.281(-2) 


4.6 


0.576(- 


-2) 


-0.198(- 


-2) 


-0.108(- 


-2) 


-0.115( 


-3) 


0.231( 


-2) 


-0.326(- 


-2) 


0.185(-2) 


5.0 


0.272(- 


-2) 


-0.929(- 


-3) 


-0.584(- 


-3) 


-0.658( 


-3) 


0.143( 


-2) 


-0.186(- 


-2) 


0.113(-2) 



current in Eq. (j2.4p ; with NLO the OPE currents in Eq. p.7|) ; with N2L0 the relativistic correction to the one-body 
current in Eq. (|2.10l) : with N3LO(loop) the one-loop current in Eq. (12. 29^ : with N3L0(min) the "minimal" contact 
current in Eq. ^T^] with N3L0(nm) the "non-minimal" contact current in Eq. ^T^] and finally with N3L0(0PE) 
the OPE currents at N3L0 given in Eqs. ^pJ2n\ and ((2?24| . The NLO and N3LO(loop) are purely isovector, while the 
remaining operators have both isoscalar and isovector terms. As in the case of the charge form factors, the expected 
suppression of the NnLO corrections as (g/A^)", where we have taken Q ^ q as the "low-momentum" scale and 
A^ = 700-800 MeV as the chiral-symmetry breaking scale, does not appear to be satisfied (not even at the smallest 
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TABLE XIV: Same as in Table IXlIIl but for ^'H. 



q(fm-i) 


LO 


NLO 


N2LO 


N3LO (L) 


N3LO (min) 


N3LO(nm) 


N3LO(OPE) 


0.2 


2.56 


0.199 


-0.325(-l) 


0.413( 


-1) 


0.178(-1) 


0.659(-l) 


0.517(-1) 


0.6 


2.11 


0.188 


-0.319(-1) 


0.380( 


-1) 


0.163(-1) 


0.612(-1) 


0.439(-l) 


1.0 


1.47 


0.164 


-0.290(-l) 


0.333( 


-1) 


0.138(-1) 


0.530(-l) 


0.312(-1) 


1.4 


0.894 


0.131 


-0.232(-l) 


0.272( 


-1) 


0.107(-1) 


0.432(-l) 


0.176(-1) 


1.8 


0.477 


0.964(-l) 


-0.159(-1) 


0.211( 


-1) 


0.763(-2) 


0.334(-l) 


0.627(-2) 


2.2 


0.221 


0.656(-l) 


-0.894(-2) 


0.156( 


-1) 


0.502(-2) 


0.247(-l) 


-0.119(-2) 


2.6 


0.834(-l) 


0.418(-1) 


-0.354(-2) 


0.111( 


-1) 


0.299(-2) 


0.176(-1) 


-0.450(-2) 


3.0 


0.188(-1) 


0.252(-l) 


0.117(-3) 


0.764( 


-2) 


0.155(-2) 


0.122(-1) 


-0.615(-2) 


3.4 


-0.591(-2) 


0.145(-1) 


0.156(-2) 


0.508( 


-2) 


0.607(-3) 


0.815(-2) 


-0.576(-2) 


3.8 


-0.117(-1) 


0.796(-2) 


0.201(-2) 


0.327( 


-2) 


0.353(-4) 


0.531(-2) 


-0.469(-2) 


4.2 


-O.lOl(-l) 


0.419(-2) 


0.177(-2) 


0.203( 


-2) 


-0.275(-3) 


0.337(-2) 


-0.350(-2) 


4.6 


-0.666(-2) 


0.210(-2) 


0.128(-2) 


0.121( 


-2) 


-0.415(-3) 


0.207(-2) 


-0.243(-2) 


5.0 


-0.365(-2) 


0.993(-3) 


0.777(-3) 


0.681( 


-3) 


-0.451(-3) 


0.124(-2) 


-0.159(-2) 




q (fm"') 

FIG. 11: Cumulative contributions to the '^He and magnetic form factors, obtained with the N3LO/N2LO Hamiltonian and 
cutoff A = 500 MeV, from the components of the current operator order by order. Set III is adopted for the isovector LEC's 
dY and dX in Table Hill 



q values). 



V. CONCLUSIONS 

In the first part of this study fSecHIland Appendices 1X1 and [Bjl . we have clarified the origin of some of the differences 
in the N3L0 and N4L0 corrections to the current and charge operators, reported in Ref. 0, Q and in Ref. In 
contrast to the authors of Ref. [3], we have not yet provided a complete derivation of the contributions associated 
with loop corrections to tree-level (OPE) current and charge operators (although some were discussed in Ref. [1|); in 
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particular, we have not carried out a full-fledged renormalization of these operators in our formalism. However, as 
pointed out in Sec. Ill Al the renormalized OPE current in Eq. (4.28) of Ref. [| leads to the same magnetic moment 
operator obtained from the currents in Eqs. (j2.7|) . ()2.22|) . and (j2.24p of the present work, with the understanding, of 
course, that the EEC's entering these equations are assumed to have been renormalized. There remain differences in 
the pion-loop corrections to the short range charge operator, Eq. (5.5) of Ref. [5] and Eqs. (|2.39|) - (|2.40|) . the latter 
presumably due to the different ways in which non-iterative pieces of reducible contributions are isolated in the two 
formalisms. The authors of Refs. [3, ^ use TOPT in combination with the unitary transformation method '63] to 
decouple, in the Hilbert space of pious and nucleons, the states consisting of nucleons only from those including, in 
addition, pious. In contrast, we construct a potential such that, when iterated in the Lippmann-Schwinger equation, 
leads to a T-matrix matching, order by order in the power counting, the xEFT amplitude calculated in TOPT j3l,[65j. 

In the second part of this study, we have provided predictions for the static properties, including charge and 
magnetic radii and magnetic moments, and elastic form factors of the deuteron and trinucleons. The wave functions 
describing these nuclei were derived from either xEFT or conventional two- and three-nucleon potentials. The matrix 
elements of the X^FT charge and current operators were evaluated in momentum-space with Monte Carlo methods. 

The xEFT calculations (based on the N3L0 potential) and the hybrid ones (based on the AV18) reproduce very 
well the observed electromagnetic structure of the deuteron for momentum transfers q up to 2-3 fm~^. In some cases, 
as in the A{q) structure function, the agreement between the experimental and ^EFT calculated values extends up to 
q < 6 fm~^, a much higher momentum transfer than one would naively expect the present expansion to be valid for. 
On the other hand, the measured B(q) structure function is significantly under-predicted (over-predicted) for q ^ S 
fm^^ in the xEFT (hybrid) calculations. The xEFT results, in contrast to the hybrid ones, have a rather large cutoff 
dependence. This cutoff dependence originates, in the hybrid calculations, solely from that in the N3L0 current, 
while in the xEFT calculation it also reflects the A dependence intrinsic to the potential (the N3L0 for A — 500 MeV 
or N3L0* for A = 600 MeV). 

The calculated ^He and ■^H charge form factors are in excellent agreement with data up to q < 3 fm^^. However, 
the observed positions of the zeros are not generally well reproduced by theory, and the measured '^He (^H) form 
factor in the region of the secondary maximum at g ~ 4 fm~^ is underestimated (overestimated) in both xEFT and 
hybrid calculations. A glance at the -Fc((?) and F^{q) in Fig. [5] suggests that two-body isovector contributions to the 
charge operator should be considerably larger (in magnitude) than presently calculated, in order to shift the zero in 
F^{q) to smaller q. 

The isovector currents at N3L0 depend on two EEC's {dX and d^), which have been fixed in one of three different 
ways: by reproducing the experimental np radiative capture cross section cr„p and isovector magnetic moment 
of the trinucleons simultaneously (set I); by using A dominance to constrain and by determining dX so as to fit 
either CT„p (set II) or (set HI). Set I is not seriously considered for the reasons explained in Sec. lIV Al The '^He and 
•^H magnetic form factors calculated with N3L0 currents corresponding to set HI, while in excellent agreement with 
data for q i fm~^, under-predict them at higher momentum transfers. On the other hand, set II N3L0 currents 
in the xEFT calculations (based on the N3LO/N2LO and N3LO*/N2LO* Hamiltonians) would lead to significantly 
better agreement with data over the whole range of momentum transfers (see Fig. 1101) , but would overestimate the 
observed fiy by ~ 3%. 
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Appendix A: Minimal contact currents 

In this appendix, we show the equivalence between the minimal contact current in Eq. (3.11) of Ref. |2j] and that 
given in Eq. (|2.20l) in terms of the known low-energy constants (EEC's) Ci, . . . , Cy. One way to achieve this is to start 
from the Lagrangian given in Eq. (2.13) of Ref. 1661] (additional terms with fixed coefficients proportional to 
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have been ignored) 

L = -\csOs ~\ctOt ~\cx{Px + 202) 
+ ^^2(2 O2 + O3) " ^^3(09 + 2 O12) 

+2 Olo)- 4^7(07 + 08 + 2 O13) , (Al) 
lb 

where the operators Oi are the standard set in Table I of Ref. [1^, and then to gauge the gradients as VA^ ^• 
VA^ — zeeAfAA^to obtain 

ji!min = ^(n.^-T2,.) (K1-K2) 
C4 

+ — - ■T2,z) o-i ■ a-2 (Ki - K2) 

i C5 

— (cTi + cr2) X (ei ki + 62 k2) 

+ -3- (n.z - T2,2) [cTi • (Ki - K2) (T2 
o 

+(T2-(Ki-K2)*Ti] . (A2) 

Because of the antisymmetry of two-nucleon states, we have 

jiln^-^^^^^^'^^^'^jiln (A3) 



where p^v^'^°^ P'^ ^ and P'^ are, respectively, the space, spin and isospin exchange operators. Making use of the 
identities: 

p^pacc(Ki_K2) =-(ki-k2)/2 , (A4) 

P"*Ti •0-2 = (3-cri •o-2)/2 , (A5) 

-P"(TM-r2..) = i(xi XX2). , (A6) 



p ^ p space 



cri-(Ki-K2) Cr2 + (T2-(Ki — K2) (Ti 



1 r 

2 



(ki-k2)(l-<Ti-<T2)+cri-(ki-k2) cr2+<T2-(ki-k2) (Ti , (A7) 



Eq. (|2:20l) in the text follows. 

An alternative way to proceed is to express the original set of LEC's C(, . . . , C[/^ entering Eq. (3.11) of Ref. 
in terms of the twelve independent LEC's Ci, . . . , C7, Cj*, . . . , C| (in the notation of Ref. 0), and then to set the 
C* = 0, that is, to ignore the currents induced by these terms, since they are suppressed by 1/m^. Substituting 

1 — ' '^2 - (--l ^ J <-'3 - ^g^2 , 



c' = --a. 



1 



2'^ii ^ + ~C'7 , Cg 



16 



1 



1 



12 



C3 



^13 



14 



(A8) 



into Eq. (3.11) of Ref. % we find 

Eq. (3.11) of Ref. [1] = Eq. (A2) - z e C4 (ei + 62) 

x(cri +0-2) X q , (A9) 

and the difference can be absorbed into a redefinition of C[c^, since (T;f +r|)(cri +0-2) x q = after antisymmetrization. 
Notice also that, in view of the identity O4 + O5 = Og + O15 (which was derived in Ref. [g^, apart from O15, of 
no relevance there), among the operators of the sub- leading contact Lagrangian, the operator Og is redundant, and 
indeed the dependence on the associated LEG Cg cancels in the observables. 
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FIG. 12: Time-ordered diagrams illustrating one ol the classes ol loop corrections to the single-nucleon current. Nucleons, 
pions and photons are denoted by solid, dashed, and wavy lines, respectively. 




Appendix B: One-loop short-range current and charge operators 



In this appendix, we discuss the contributions associated with panels (h)-(k) in Fig.[2]for the current operator, and 
(g)"(j) in Fig. [3] for the charge operator. We begin with the current operator. The contributions of diagrams (h) and 
(j) in Fig. [2] vanish, while the contribution of diagrams of type (i) was obtained as (conventions for q-intcgrations and 
(5-functions are the same as in Ref. 

.(1) ^ 2z ^4?^ ^ ^2). / %i + q2 - q) 

xD(cji, W2) (qi - q2) cTi • q2 0-2 • qi , (Bl) 



where 



' UJ1UJ2 + w| 



UJ1UJ2 (UJl + i02) 



Before analyzing diagram (k), we need to consider the leading and next-to- leading contributions to the single-nucleon 
diagrams shown in Fig. 1121 For simplicity, we define the vertices 

Vi = cr • qi Tq ,¥2= a ■ci2n , (B3) 



Vj = -ieeabz (qi - q2) • 



(B4) 



Then the current reads 



ViV2V-t 

AuJlLU2 



1 



1 



1 



Ei — LOi — E Et - uji - ijJ2 — E 

1 



Ei — LOi — E Ei — Lo^ — UJ2 — E 

1 1 



E^ - UJ^ - UJi - UJ2 - E' Ei - UJ-, - iL!2 ~ E 



(B5) 



where Ei = E + lj^ , and E is the energy of the intermediate nucleon of momentum p. After expanding the energy 
denominators as in Eq. (|2.18p to include linear terms in the nucleon kinetic energies, we find, up to next-to-leading 
order included, 



ViV2V^ 



UJtUJ. 



+ D'{uJi,uj2){E' -E) 



r"^2 



-|-I?(a;i,a;2)(£' - E) + D^{wi,i02) uj^ 



(B6) 
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FIG. 13: Set of time-ordered diagrams for tiie contribution illustrated by the single diagram (k) in Fig. [2] Notation as in 
FiE.fll 



class A 



El 



class B 



class C 



where 



and 



D(a;i,W2) 
D^(a;i,W2) 



UJi +2UJ2 

wj'wKcji + UJ2) 

UJl ~ LO2 

UJ1UJ2 (wi + 0J2Y 



D^D' with qi ^ q2 , D + D' ^ 2 D 



(B7) 
(B8) 
(B9) 

(BIO) 



We now proceed to analyze the contributions of diagrams of type (k) in Fig. [2] To this end, we show in Fig. [13] 
the complete set of time-ordered diagrams of the same topology as (k) , which we have separated for convenience into 
the three classes A, B, and C. Class A consists only of irreducible diagrams, which at order n = 1 or eQ, i.e., in the 
static limit, lead to 

class A = - ^ KyViVcTV2i5(a;i,w2) 



-VjViVctV2[D{uji,uj2) + D'{uji,uj2)] , 



(Bll) 



where the vertices Vi, V2, and V-y are defined as above (with the spin and isospin matrices now referring to nucleon 
1), and 



VcT — Cs + Ct (Ti ■ (72 

On the other hand, to order e Q included, class B gives 



(B12) 



class B = 



1 



4wiW2 E,-E^-E2 
E\ — E\ \ 1 



1 



UJi{lJI + UJ2) 



Ei+ijj^- El 



UJl 



0J2 



UJ2{ljJl + W2) 



Ei — El — E2 

UJl + UJ2 



UJl + UJ2 

El — El 



1 



1 



Ei+uj^- El 



UJl -I- UJ2 



UJ2 



UJ1UJ2 \ UJl 

_ V-,VctViV2 1 

A:UJlUJ2 UJ2{UJI + UJ2) 



(B13) 
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FIG. 14: One-loop correction to the OPE current (only one among the possible time-orderings is shown). Notation as in Fig. 1121 




where we have used energy conservation between the initial and final states Ei + E2 + uj-y — E[ + E2. We now note 
that the irreducible contribution from the last diagram (in class B) is cancelled by the second term in the next to last 
line of the above equation, so that we are left with 



class -^-'^-^^^^ 



4 E,-E[-E.. 



2 

2-3 + ujj D^{uji,uj2) + {El - El) D{uji,uj2) + {El + cj^ - Ei) D' {uji,uj2) 



l'^2 



1 KV(-tV^iV9 

= VcT =^ h.. + ^ , ^ K>^2) , (B14) 

E,-E^-E2 4 

since Ei+uj^ — Ei = {Ei — Ei — E2) + {E^ — Ei) and V-^ commutes with each of the remaining vertices. The first term 
represents an iteration, while the recoil-corrected class B contribution is simply given by VyVcTViV2 D' {uji,ijJ2)/4:. 
A similar analysis for class C leads to the recoil-corrected class C contribution given by V^ViV2Vct D{uj 1,102)/ 4:. 
Therefore combining the contributions from classes A, B, and C, we find 

i'^ - \v^Vi [V2,Vct] D{ui,L02) 

+ \ [VcT , Vi] V2 D'{lui,lo2) , (B15) 

or explicitly 

j(i) ^ 2zi^|^ri, / 5{tii+ci2-q)D{uJi,uJ2) 

x(qi - q2)cr2 • q2 X qi + 1 ^ 2 . (B16) 

The currents jp'' and j^'' obtained here are in agreement with those in Eq. (5.2) of @, but for an overall factor of 2. 
Ultimately, this difference has no impact, since 

jP^ + Jk^ oc (ri X T2)z (cTi X 0-2) X q 

+ (2Ti,(T2 + 2T2.cri) xq=0 , (B17) 

which vanishes due to antisymmetry of the two-nucleon states. 

In Ref. 01 we had not considered the next-to-leading order contributions to the single- nucleon jmr vertex when 
deriving the one- loop correction to the OPE current shown in Fig. [T3] (see Appendix E of Ref. [2] ) . As a consequence 
we had failed to isolate the correct non-iterative piece, which had led, in particular, to a non-hermitian operator. We 
find that this term is now given by 



JOPE,loop 



X T2) 



(T2 ■ k2 



■ (Ti • k2 (g^ - 



UJ^U}^{U}+ + UJ-) 

1^2 , 



(^CTi • q k2 • p CTi • p k2 • q 



(B18) 



where w± — y^q±pp + 4TO^. 

Next, we turn our attention to the charge operator. In Ref. Q we showed that the contributions of diagrams (g)-(h) 
in Fig. [3] vanish. However, in light of the previous considerations, those due to diagrams (i) and (j) given there need 
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FIG. 15: Set of time-ordered diagrams for the contribution illustrated by the single diagram (i) in Fig.|3l Notation as in Fig. 1121 
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'7^ 'X 



(d) 



to be revised. Indeed, an analysis similar to that carried out above leads to the single-nucleon charge operator (see 
Fig. [T^ up to next-to-leading order included 



ViV2V^ 



+ D'{uJi,uj2){E' -E) 



+D{LUl,LU2)iE -E)+ D^{LUl,UJ2)i^'^ 



and to a contribution for diagram (j) (see Fig. [131) which reads 



We have defined K, = — z e eabz , and 
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(B19) 
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We now revise the derivation of the charge operator of type (i) illustrated in Fig. [3] The associated time ordered 
diagrams are represented in Fig. I15[ and have been separated into three classes. In Ref. [Sj, the expression reported 
in Eq. (53) has been obtained by accounting for the recoil corrected class A diagrams only, i.e. 
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FIG. 16: Diagrams for the single-nucleon contributions. Notation as in Fig. 1121 
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where Vi and Vct are given in Eqs. (|B3j) and ()B12p . respectively, and V{ and are defined as 



(B26) 



Classes B and C involve only reducible diagrams. First, we analyze the single-nucleon contributions entering the 
reducible diagrams, represented in Fig. 1161 We account for leading, next-to- leading, and next-to-next-to-leading order 
corrections in the expansion of the energy denominators. The contributions Af^r and p-y, associated with panels (a) 
and (b), respectively, read 
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(B27) 
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where the energies are as indicated in Fig. 1161 and energy conservation (E + uj-y — E') has been used. In terms of 
these, the contributions of diagram (a) and (b) in class B are given by 
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(B30) 



where we have identified and isolated the nucleon self-energy terms (to be reabsorbed by mass counter-terms). Eval- 
uation of panel (c) leads to 
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(B31) 



where the first term represents an iteration with and the contact interaction, and the second term is the recoil 
correction contributing to the two-nucleon charge operator. Finally, the contribution of panel (d), in which, in contrast 
to panels (a) and (b), the self-energy insertion is between the photon absorption and contact interaction, is expressed 
as 
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, (B32) 
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where the last term is a two-nucleon term, and we interpret p-y as a further correction to the single-nucleon 7A'' vertex. 
However, the analysis and proper interpretation of this type of corrections are beyond the scope of the present work. 
After a similar analysis of the class C diagrams is carried out, we find that the complete B+C contribution reads 

2^ \ViV:^ ViVcT - V:^Vi ViVcT 

+■^^0x^1!^^' Vi - VctViVIV^ 
J9\ 



classes B + C 



^ti,.(Cs + Ct<ti.<t2) / % 



which combined with Eq. (jB25|) leads to the type (i) charge operator in Fig. |3] 
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Appendix C: Loop integrations 



In this appendix, we outline the derivation of the two-body charge operators at one loop, listed in Sec. IIIBI For 
the sake of illustration, we consider the contribution of panel (f ) in Fig. |3l given by (in the notation of Ref. 
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which can conveniently be written as 
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- (n X T2)^ (q/2 + p) • (p - k) <T2 • (q/2 - p) x (p - k) , (C3) 

and the momentum k defined as in Eq. (j2.2l) . We now use standard techniques [67|] to express the product of energy 
denominators in the following way 



which, in terms of 
simply reads 
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where 

X'^{zi,Z2) = {zi + Z2) qV4 - [ (^1 - Z2) q/2 - (1 - zi 

-Z2)k]' + (l-zi-Z2)k2+m2 . (C7) 

After these nianipulations, the charge operator can finally be written as 

.1/2 

/ nnf nr I 

[p'^ + \\x,y)]-'' + l^2, (C8) 



^ f\xx f'\v f iV'(q,k,p') 

TT Jq J-1/2 Jp' 



X 



where the function N' is obtained from N by expressing p in terms of p' via Eq. (jC5|) . We have also changed variables 
in the parametric integrals by introducing [67 1 

x^zi + Z2, xy ^ {zi - Z2)/2 , (C9) 

such that 

rl-zi 1-1/2 

dzi dz2 — > dxx dy . (CIO) 

Jo Jo J-1/2 

The function N' is a polynomial in p', and the p '-integrations are carried out in dimensional regularization (see 
App. A of Ref. [I])- They are finite and lead to the charge operator given in Eq. (|2.38p . 
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